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Abstract 

A theoretical model for the calculation of the phonon-drag thermopower, S^, in degenerately 
doped semiconducting single-wall carbon nanotubes (SWCNTs) is proposed. Detailed calculations 
of are performed as a function of temperature, tube radius and position of the Fermi level. 
We derive a simple analytical expression for that can be utilized to determine the free carrier 
density in doped nanotubes. At low temperatures shows an activated behavior characteristic of 
the one-dimensional (ID) character of carriers. Screening effects are taken into account and it is 
found that they dramatically reduce the magnitude of . Our results are compared with previous 
published experimental data in bulk p-doped SWCNT materials. Excellent agreement is obtained 
in the temperature range 10-200 K for a consistent set of parameters. This is a striking result in 
view of the complexity of these systems. 
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I. INTRODUCTION 



Thermopower, 5", is an important transport coefficient that offers valuable information 
about the electronic structure, the scattering processes and the mechanisms of carrier- 
phonon coupling in a system. In the last few years there has been growing experimental 
interest in S of single wall carbon nanotubes (SWCNTs). Several groups have reported 
thermopower measurements on bulk SWCNT materials (e.g., mats, fibers, ffims) [liJgi and 



on individual SWCNTs 



10 



12| . However, only modest progress has been made up to now in 



understanding the unique features of S in these systems. Interesting issues concerning the 
large positive thermopower (~ 80 uV/K) in pristine samples {2, 4-6, 8|, the change of sign 
of S upon exposure to oxygen jl, 6| and the effect of carrier-phonon coupling [t-Q, 13-15| 
on S still remain open. 

S consists of two additive contributions which are diffusion, S'^, and phonon-drag, S^. 
S"^ is due to the carrier diffusion in the presence of a temperature gradient and for degen- 
erate systems varies linearly with T according to Mott's expression. originates from the 
interchange of momentum between acoustic phonons and carriers via the carrier-phonon 
interaction. The first theoretical models for the study of the phonon drag in metals [l^ 



and semiconductors 



171 ] were developed half a century ago. More recently, extensive theo- 



retical and experimental work has been carried out on of low-dimensional semiconductor 
structures 



18 



20|. 



Recent experiments on 5* in p-doped SWCNT films and fibers js, 9| provided clear evi- 
dence for the presence of S*^ at T > 15 — 20 K. On the theory level, however, there is still 



15 



So far, the 



, [l3| . However, 



an ongoing discussion about the role of in measured thermopower 
theoretical studies of 5*^ are confined to metallic armchair (10,10) tubes 
in perfect metallic tubes with mirror electron-hole symmetry both S'^ [7] and 
expected to be negligibly small compared to the experimental data, due to the competition 
between the opposite contributions of electrons and holes. We note that the accuracy of the 



ii5i 



are 



existing theoretical models 



Mahan 



for in metalhc tubes has been questioned recently by 
15| . Also, a recent theoretical work 2l| pointed out that thermopower vanishes in 



one- dimensional conductors with a linear energy dispersion (as in the case of metallic tubes) 
due to electron-hole symmetry. 

In this paper we propose a theoretical model for the phonon-drag thermopower in semi- 



2 



conducting SWCNTs that are characterized by a non-hnear energy dis per sion. (A brief 
discussion on the behavior of 5*^ in this kind of nanotubes appears in [20!].) We suggest 
that the measured thermopower in doped samples is due to the contribution of degener- 
ate semiconducting nanotubes. In our model originates from carrier-phonon intraband 
scattering within the first ID subband. As we discuss below, the dominant contribution to 
is made by long-wavelength acoustic phonons that backscatter carriers across the Fermi 
surface. In this case the carrier-phonon coupling is much weaker in metallic tubes than in 



22| and, consequently, is expected to be substantially larger in the 



semiconducting tubes 
latter ones. 

We note that upon chemical or electrostatic doping the Fermi level can be pushed into 
the conduction or valence band and the degenerate semiconducting tubes can be considered 
as one-dimensional metals. Therefore the terms "metallic" and "semiconducting" refer only 
to the different electronic structure in the two types of tubes (see, for example, Ref. |23l)]. 

There are two equivalent theoretical approaches to the problem of phonon drag 20| . In 
the first approach phonons are perturbed in the presence of a weak temperature gradient 
VT. Non-equilibrium phonons transfer part of their momentum to carriers due to the 
carrier-phonon coupling. Then the phonon-drag contribution to the thermoelectric current 
J'J = L^VT is calculated by solving the coupled Boltzmann equations for carriers and 



acoustic phonons 



13 



16 



24 



25| . The phonon-drag thermopower is readily obtained by 



5*^ = — L^/cr where a is the carrier conductivity. In the second approach carriers are 
accelerated isothermally in the presence of a weak electric field E and impart some of their 
momentum to phonons due to the carrier-phonon coupling. Then the resulting phonon heat 
current and the phonon-drag contribution to the Peltier coefficient is calculated 



17 



26- 



3l|. This method of evaluating is referred as Il-approach [17] because it provides a 



direct estimation of the Peltier coefficient. The equivalence of the above two approaches 
is secured by Onsager's symmetry relation. In this paper we follow the second approach 



which is more gen era. 



semiclassically |28l431|. 



and it can be applied even in systems where carriers do not behave 



The paper is organized as follows. In Sec. II we introduce the theoretical model for the 
calculation of S"^ in the semiclassical transport regime. An explicit expression for is 
derived in Sec.IIB and in Sec.IIC we derive a simple approximate expression for 5*^ for the 
case of a highly degenerate semiconducting tube. Numerical results for as a function of 



3 



temperature, tube radius and position of Fermi level are presented in Sec. III. In the same 
Section we discuss the effect of screening. In Sec. IV we compare our theory with available 
experimental data for acid-doped bulk SWCNT samples. 



II. THEORY 



A. Description of the physical system 



We assume that the nanotube is a long indefinitely thin cylinder of radius R and length 



L. The nanotube axis is along the 2— direction. The carrier wave function is 



32| 



jie 



6{r-R), 



(1) 



where, r is the space vector, k is the carrier wave vector along the axial direction, 6 is the 
azimuthal angle and / labels ID orbital subbands associated with the carrier confinement 
along the circumference. We assume, that the Fermi level, Ep, is located between the first 
and the second ID subbands (i.e., only the ground subband is occupied). Then, the carrier 
energy is 

Ek = Ei + ^ (2) 



2m* 

where m* is the carrier effective mass and Ei denotes the position of the first van Hove 
singularity. 

In carbon nanotubes phonons also exhibit ID character. The lattice displacement at a 



point r is 



22| 



u r 



- zqz imt 



(3) 



where, f]mq is the polarization vector, q is the phonon wave vector in the axial direction and 
m = 0, ±1, ±2, ... denotes the phonon modes associated with phonon confinement along the 
circumference. Due to the conservation of angular momentum only the three low-energy 
acoustic modes with m = (the so-called twisting, stretching and breathing modes) con- 
tribute to the carrier-phonon intraband scattering. The phonon frequencies and polarization 
vectors have been calculated within the continuum model proposed by Suzuura and Ando 



22|. 



The carrier-phonon interaction in carbon nanotubes has been studied in several texts 



within the tight-binding approximation 



33 



- l37j or a continuous elastic theory 



22 



38 



4q1- 



Here we follow the continuous model of Suzuura and Ando 



22] according to which the 



carrier-phonon coupling is described via the acoustic deformation potential 

f/(r) = D 



1 due duz Ur 

H + — 

R dd dz R 



(4) 



where D is the deformation potential constant. The deformation potential approximation 
provides a good description of the carrier interaction with long-wavelength acoustic phonons. 
The last term in Eq. (jlj) accounts for the nonzero curvature of the nanotube 22|. The 
twisting mode does not participate to carrier-phonon scattering via the deformation potential 
coupling. Moreover, in the long- wavelength limit {qR <^ 1), which is the regime of our 
interest, the breathing mode is dispersionless and does not contribute to . Thus, in what 
follows we consider only the stretching mode which is characterized by a linear dispersion 
= "^skl where, Vg is the sound velocity. The phonon polarization vector, r) = {riQ,riz,rir), 



for this mode in the limit qR ^ 1 is 



Vq = (0: 



1 —iuqR, 



(5) 



where a = a/1 + z/^g^i?^ and u is Poisson's ratio. Ignoring the terms proportional to q^R 



the above expression becomes identical with the one derived by De Martino et al. 



4l|. 



B. An explicit expression for the phonon-drag thermopower 

We assume a small electric field E in the axial direction of the nanotube. The presence 
of E creates a net flux of carriers along the axis of the tube which results in a momentum 
transfer to phonons through the carrier-phonon coupling. We calculate the resulting phonon 
heat flux Q and obtain the phonon-drag contribution to the transport coefficient 

= Q/E. (6) 

To get we utilize the Onsager's relation 

S' = ^ (7) 

where a is the carrier conductivity and T the absolute temperature. 
The phonon heat flux is given by 

Q = }^Y.h^^v,Nl (8) 
1 



where Vq = Vsq/\q\ is the phonon group velocity and A^^^ = Nq — is the first order 
perturbation of the phonon distribution function. 

The perturbation A^^^ is determined by the steady-state Boltzmann equation for phonons 
in the relaxation time approximation when VT = 0. Namely, 

where Tp^ is the phonon relaxation time associated with phonon-phonon collisions and 
phonon scattering by imperfections. For simplicity we have ignored the dependence of 
on q. {dNq/dt)ph-c is the rate of change of the phonon distribution function Ng due to 
phonon scattering by carriers. It is written in the standard form 

ph-c 1.^1^, 

-fk{l-fk')P^{k,k'), (10) 



dt 



where Qs and are the spin and the valley degeneracies, respectively, fk is the carrier 
distribution function and Pq^'^\k, k') are the transition rates at which the carrier in a state 
k is promoted to a state k' by absorbing (emitting) one phonon with wave vector q. 

When the external field E is weak Eq. (|TOl) is linearized and is solved in terms of A^"^. 
Then we get 

dNq\ _ N^q 



9t )ph~c V(?) 

gsQv sr^Y ( 

(11) 

where, Tpc{q) is the phonon relaxation time associated with scattering by carriers given by 

Tp-'iq) = 9s9.J2^k',k/[N^iN^ + 1)], (12) 

k,k' 

and Ffc/ is the average equilibrium rate of absorption of phonons with wave vector q. It is 
given by 

rk',k = mi-f,-)PfiKk% (13) 

where = f{Ek) = {exp[^{Ek-EF)] + l}-^ (with /3 = l/fc^T) is the Fermi-Dirac function 
and Pq^{k, k') denotes the transition rate in equilibrium. 



Assuming that phonon-phonon scattering and phonon scattering by impurities dominate 
over the phonon- carrier scattering (rpc ^ Tph), Eqs.Q and f|TT]) give 

,,^1 _ 9s9vTph \^ -p f fk fk' \ 

In the above equation fl is the first order perturbation of the carrier distribution function. 
It is worth noting that Eq. f|T4|) can be regarded as a starting point for the calculation of 
in all the problems treated within the Il-approach 20|. Now, by substituting the phonon 

perturbation into (|8]) we take for the heat fiux 

«-ltE-.^.r.,(-J^--J^). (15) 

To determine the perturbation of the carrier distribution function fl entering Eq. ( |T5|l 
we use the ID steady-state Boltzmann equation 

where e is the carrier charge and the RHS of Eq. (|T6l) is the rate of change of the carrier 
distribution function due to elastic collisions with static imperfections. In the relaxation 
time approximation this term is written as —fl/T{Ek) where T{Ek) is the carrier relaxation 
time. Equation (fT6|) is linearized to give 

fl = -eET{E,)v, (^^) (17) 

where, Vk = {l/h)'V]iEk = hk/m* is the carrier group velocity. 

By substituting Eq. (fT7|l into ( TT^ and making use of Eqs.(jn]), (El) and (IT^ we finally get 



S9 



aLkeT^ 



huJqVq[T{Ek)Vk - t{E, 



k'jVk' 



k,k' ,g 



xfki^-fk')Pfik,k'). 



(18) 



The above expression is equivalent to the expression derived by Kubakaddi and Butcher 25] 



25] followed a different 
16| and they calculated 



for a quantum wire coupled to 3D phonons. The authors in Ref. 
approach than this described here. They followed Bailyn's theory 
the phonon-drag contribution to the thermoelectric current that originates from the carrier 
scattering with non-equilibrium phonons in the presence of a small temperature gradient 



across the wire. Their calculation was based on the solution of the coupled equations for 
electrons and phonons. 

The transition rate Pg^{k,k') is calculated by using Fermi's golden rule. The lattice 
displacement for the stretching mode is written in second quantized form 



"(r) = E y ^ i^'^^'"'^^ + %^"''<) ' (19) 

where and are the phonon creation and annihilation operators, respectively, A = 2ttRL 
is the nanotube surface area and p is the mass density. For the stationary carrier states 
considered here one easily finds 

P^'ik, k') = ^^' J^fT^ iE,, -E,- hw,) 4',fc+, (20) 

where, = [exp{f3hu!q) — 1]^^ is the phonon distribution in equilibrium, \Uq\'^ is the square 
of the carrier-phonon matrix element for the deformation potential coupling and e(|g|,T) is 
the ID static dielectric function. By utilizing Eqs. (jl]) and ([5]) the matrix elements |f/gP in 
the limit qR <^ 1 are written as 

2 ^ (21) 

where S = D(l — z/). We note that the g-dependence of |f/qP is typical for the carrier 



interaction with longitudinal acoustic phonons via an isotropic deformation potential 
A similar expression to the one we derive here is given in Ref. 



42|. 



The dielectric function for a ID gas confined to the surface of the carbon nanotube is 



calculated within the random phase approximation 32|, |43| . For the carrier wave functions 
considered here we obtain 

e(|g|,r) = i + %^^iro(|g|i?)/o(|g|i?)M(|g|,r) (22) 

n^TTC;, \q\ 

where Jq and Kq are the modified Bessel functions of the first and the second kind, respec- 
tively, and efe is the background dielectric constant. M{\q\,T) is the standard factor that 
accounts for finite temperature effects on the static polarization function 43 



To obtain an explicit expression for we substitute Eq. (!20|) into (fT8|) . Then the summa- 
tion over k' is readily carried out by replacing /c' by /c + g as a consequence of the momentum 

8 



conservation condition imposed by the Kronecker symbol Sk',k+q- Moreover, the summations 
over q and k are transformed to the integrals 

The presence of the (5-function in Eq. ( 120|) allows the immediate evaluation of the k- 
integration. We see by inspection that 

2m* 

6{Ek+q -Ek- hw,) = ^^^C^k + qTqo) (24) 

where, go = 2vsm*/h. The minus and the plus signs correspond, respectively, to positive 
and negative q. 

Now, after some algebra, we finally obtain 

= rd<l-,7jh^.7&N^m (25) 



2-nepRkBT'^ Jo ^ e^{\q\,T)2kF " 

where, Iph = VgTph is the phonon-mean-free path, kp = \lm*{EF — Ei)/fi?Y^'^ is the Fermi 
wave number and /(g) is the product of the Fermi occupation factors 

I{q) = f{E,)[l-f{E, + nw,)] (26) 

with = |go — g|/2. In deriving Eq. ( 125|) we have ignored the energy dependence of the 
carrier relaxation time and in Eq. ( ITSll we have replaced rf-E/j^) by its value at the Fermi 



level, Tp. This is a good approximation when tvjjq <^ Ep [45|. Moreover, we have replaced 
a by ne^Tp/m* where n = gsQukp/TC is the density of carriers per unit length. Interestingly, 
becomes independent of the carrier relaxation time. 

C. An approximate expression for 5^ 

At low T and assuming that fvjjq is a small quantity compared to Ep the product /(g) is 
approximated by [24 1 

/(g) ^ nujq{Nl + l)6{Ek - Ep), (27) 
with k = \qo — q\/2. The (5-function can be written in the following form 

2m* 

5{Eu - Ep) = ^[Siq - go - 2kp) + 5{q - go + 2kp)]. (28) 



9 



We see that 5{Ek — Ep) resonates at g = go + for positive q. When the expression (|2j 
for I{q) is substituted into Eq. fl25|) the integration over q is carried out straightforwardly 
by using the condition g = go + 2fc_F- We note that go ^ 2kp and consequently, stretching 
phonons with q = 2kp make the dominant contribution to . 

Equation f l2^ is now significantly simplified and is written in the convenient approximate 
form 

S9 = 1^ ^ r (oq) 

e^{2kF,T) {e^^^^F -1)2 ^ 

where C is given by 

^ ^ 2{m*yE%hUJ2kp .^Q. 
nhepRkB 

In the above equations, U2kp = Vs2kp is the frequency of a stretching phonon with g = 2kp 
and e{2kp,T) is an approximate expression for the dielectric function. To obtain e{2kp,T) 
we replace g by 2kF in the denominator and in the arguments of the modified Bessel functions 
Jo and Kq in Eq. f l22|) . The factor M(|g|,T) is replaced by the average M{2kp,T) that is 
given by the expression 

M{2kp,T) has been evaluated numerically for several values of kp and R and we find that 
in the degenerate limit and when Tp > 5fiJUJ2kp/kB (where Tp = [Ep — Ei)/kB is the Fermi 
temperature) the following expression provides a very good fit 

M{2kp, T) = In (^^t£±3L^ [ai - as exp(-a3a;)] (32) 

where, x = Phuj2kp, «i = 1-175 ± 0.002, = 0.60 ± 0.01 and = 0.41 ± 0.01. At low 
T the effect of screening is severe and unity can be neglected in Eq. (!22|) . In this case the 
T-dependence of the dielectric function is described by Eq.( l32i) . 

At temperatures where (3fvjj2kp ^ 1 the dielectric function shows a weak T-dependence. 
Then follows the law 

oc ^e-^^'^^''F . (33) 
This activated behavior is characteristic in ID systems where the Fermi surface consists of 



two discrete points ±kp |46l. |47|. 



We note that when screening is ignored and the phonon mean-free path is constant the 
T-dependence of given by Eq. (!29|l is similar to what predicted by Scarola and Mahan 

10 




T(K) 



FIG. 1: against temperature for a SWCNT of radius 0.5 nm. Results are shown for four values 
of Ep — Ei: 60 meV(dotted line), 90 meV (dashed-dotted line), 120 meV (dashed line) and 150 meV 
(solid line) when screening is taken into account (a) and when screening is ignored (b). The phonon 
mean-free path is taken to be 1 /um. The inset shows the ratio A = / Sappr as a function of TpjT . 



|13l | for an armchair (10,10) metallic SWCNT due to interband electron scattering between 
the two linear bands. However, the absolute magnitude of in a metallic tube is expected 



to be much lower than that predicted in Ref.[13[ due to the competing contributions of 
electrons and holes to the thermoelectric current. 



III. NUMERICAL RESULTS 



We assume that the free carriers are holes and we examine the dependence of on 
temperature, the radius of the nanotube and the position of the Fermi level with respect to 
the position of the first van Hove singularity. The analysis is the same for the case of electrons 



with the only difference being the sign of . The va^ 



9v 



2, D = 24 eV 



22 



ues for the material parameters used 
48|, u = 0.2 [491, efe/47reo = 2.4 32|, 



in the calculations are Qs 
p = 3.8 X 10^^ Kgr/m^ and = 19.9 km/s 22|. The hole effective mass is taken to be 
m* = 'me/22.7 R where R is the tube radius in nm 49|. We assume that Iph = 1 pm. 



11 



0,00 0,(32 0,(34 0,06 0,08 0,10 0,12 

FIG. 2: Dielectric function at q = Ikp as a function of the inverse temperature for a SWCNT of 
radius 0.5 nm. The soHd, dashed, dashed-dotted and dotted Hnes correspond to kp = 0.4, 0.45, 
0.5 and 0.55 nm~^, respectively. 

In Fig. la we show the evaluated from Eq. fl25l) for a p-type SWCNT of radius 0.5 nm as 
a function of T. The sohd, dashed, dashed-dotted and dotted hnes correspond, respectively, 
to Ep — El = 150, 120, 90 and 60 meV. We note that at temperatures where the carriers 
are non-degenerate we have taken into account the thermal broadening effects on cr. To 
assure the accuracy of the approximate expression (12^ . in the inset of Fig. 1 we show 
the ratio A = / S^^^^ as a function of Tp/T for Ep — -Ei = 60 meV. and S'f^^^ are 
calculated from Eqs. fl25l) and fl29l) . respectively. Calculations of A for 90, 120 and 150 meV 
also fall on to the same curve. We can see that in the degenerate limit the approximate 
result agrees very well with the exact expression for . Finally, in Fig. lb is calculated 
in the absence of screening, e(|g|,T) = 1. It turns out that screening induces a strong 
suppression of by 1-2 orders of magnitude. Inspection of Eq. fl2^ shows that screening 
effects become more severe as R decreases. We note that in the absence of screenin g 



levels off at high T in agreement with previous estimations in metallic SWCNTs j8|, Il3 |. 
However, when screening is introduced 5*^ shows a quasi-linear T-dependence at high T due 
to the temperature dependence of the dielectric function. The dielectric function e{2kF,T) 
as a function of the inverse temperature for a SWCNT with R = 0.5 nm is shown in Fig.2 

In Fig. 3 we show the dependence of 5*^ on the Fermi level with respect to the position of 
the first van Hove singularity for temperatures 50 < T < 300 K. The shown structure is due 
to two competing mechanisms which are the suppression of the carrier-phonon scattering 
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FIG. 3: as a function of Ep — Ei for various temperatures. The phonon mean-free path is Ifim 
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FIG. 4: as a function of the nanotube radius for various temperatures. Iph = 1 /im. At 
T > 100 K fohows approximately a R~^-^ law. 

and the increase of l/e(g) as kp increases. The tube radius is 0.5 nm. 

Finally, in Fig. 4 we present the calculated values of as a function of the nanotube 
radius. At temperatures higher than 100 K we find that 5*^ follows a law close to oc R~^'^. 
At lower temperatures shows a weaker dependence on R especially at large values of R. 



IV. COMPARISON WITH THE EXPERIMENT-DISCUSSION 

So far there is no clear evidence about the phonon-drag effect in isolated SWCNTs. The 
most relevant experiments were performed by Yu et al. |lli] in an individual SWCNT at 



13 



temperatures above 100 K. The observed thermopower showed a hnear T-dependence which 
was attributed to the hnear diffusion component and a constant phonon-drag component of 
about 6 ^V/K without, however, excluding the possibihty of an additional contact effect. 
According to our analysis in Section III, the phonon-drag thermopower at relatively high 
temperatures shows a quasi-linear T-dependence and this makes difficult the separation of 
the diffusion and the phonon-drag contributions. Nevertheless, when the calculated values 
for shown in Fig. la are fitted by a linear function of T we find that the intercepts vary 
from 1.4 to 7.4 /iV/K when the position of the Fermi level with respect to the first van Hove 
singularity varies from 60 to 150 meV. These values are in agreement with the experimental 



estimation of in Ref. 



ll| . We note that the intercepts depend linearly on the phonon 

lati 




mean-free path and vary approximately as R~^'^. 

Vavro et al. js] and Zhou et al. [q] have reported thermopower measurements in p-doped 
bulk SWCNT samples in a wide temperature range (10-200 K) that show clearly the sig- 
nature of phonon drag. Normally, in bulk samples nanotubes are self-organized into long 
"ropes" , which contain a large number of nanotubes (tens to hundreds) [50| , forming a 3D 
network of complex geometry. Thermopower in these nanotube networks exhibits a very 
similar behavior as this of an individual nanotube described in Section III. We have recently 
proposed a simple argument based on a model of parallel conductors which suggests that in 
a network with homogeneous doping and with a narrow distribution of tube diameters the 



measured thermopower resembles that of an individual tube 



20y. The resistivity measure- 



ments in the samples under consideration showed weak coupling between metallic nanotubes 



Is 



9| and hence the contribution from metallic tubes to the total conductivity is neglected. We 
also recall that the contribution of metallic tubes in S is expected to be small compared to 
this of semiconducting tubes. Therefore, we can use the theory for isolate semiconducting 



SWCNTs developed here to interpret the data in 



In Fig. 5 the circles are the measured thermopower for a bulk sample prepared by pulsed 
laser vaporization (PLV) and doped with HNOajol. The tube radius is i? = 0.68 ± 0.04 nm. 
At low temperatures (up to 100 K) we fit the data for the total thermopower, by the 
expression 



S=—^— — ^— 5T + ATil-B\nT). (34) 



C 1 e^'^^fc 
T^e\2kF,T) (e/5'""^2V - 1) 
The first term is the approximate expression (129|) for and the second term corresponds 

to the diffusion component S'^. The sample is highly degenerate and at temperatures up 
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to 100 K Eq. (l29|) accurately describes S^. The T-dependence introduced by the dielectric 
function is given by Eq. f l32|) . The values we obtained for the parameters kp, A and B are 
shown in Table I. 

The logarithmic term in S"^ secures an excellent fit to the measured thermopower at 
all temperatures up to 100 K. If this term is neglected the theoretical values for the total 
thermopower are significantly larger than the experimental ones at high temperatures. We 



speculate that the TlnT term in S'^ is due to 2D weak localization (WL) effects 51 |. 
If this speculation is valid we would also expect a signature of WL in the conductivity 
measurements. We note that the relative change in conductivity should be the same as in 



S'^ but with an opposite sign [51|, l52[. Interestingly, we find that at temperatures 10 < T < 



100 K the conductivity follows the law 

a = (7o(l + 5'lnT), (35) 

where the value of B' is shown in Table I. We see that B and B' agree to each other 
remarkably well. The origin of the 2D WL in these samples is not well understood. It is 
likely related to the individual rope although in this case a ID localization behavior would 
be expected [53i]. However, the phase coherence length, in the samples we discuss here, is 



coniparable to the diameter of the rope [5j] and the 2D limit might be approached. Langer et 
al. |55[ have also observed a In T dependence of the conductance for an individual multiwall 
CNT at 0.1-100 K which was attributed to_2D WL. Finally, we should remark that WL is 
expected to have a negligible effect on [5g]- 

Now, by using the values for kp we obtained from the fitting of the thermopower data 
at low T we calculate in the whole temperature range from 10 to 200 K by using the 
exact expression fl2^ . The only remaining unknown is the value of the phonon- mean- free 
path Iph which is determined from the experimental data when the diffusion contribution is 
subtracted. We find that Iph = 0.6 nm. This value is consistent with the values 0.25-0.75 fim 



reported recently for an individual SWCNT 



ll| . Our estimation for the total thermopower 



is shown as solid line in Fig. 5. The dashed and the dashed-dotted lines correspond to the 
phonon-drag and the diffusion contributions, respectively. 

By following a similar procedure as this described above we have interpreted the ther- 
mopower data for another bulk sample prepared by high pressure decomposition of CO 
(HiPco) and doped with H2SO4 js]. The tube radius varied from 0.4-0.7 nm. Conductivity 
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FIG. 5: Thermopower versus temperature. The circles are the experimental data for a bulk 
sample prepared by pulsed laser vaporization and doped with HNOajo]. The solid line is the total 
thermopower of an individual nanotube obtained as explained in text. The dashed and the dashed- 
dotted lines correspond to the phonon-dras and the diffusion contributions, respectively. In the 
inset the circles are the conductivity data [sl and the solid line is the fit by using Eq. ([35]) . 



measurements for this sample (designated as HPR93C) appear in [57|. The experimental 
data for the ratio S/T are shown as squares in Fig.6. The values for the fitting parameters 
kp, A and B are shown in Table I. We also present the value of B' for comparison. In the 
calculations the tube radius is taken to be the average i? = 0.55 nm while for the phonon- 
mean-free path we obtained the value Iph = 0.4 nm. The calculated values for S/T is shown 
as solid line in Fig.6. 

TABLE I: The values for the parameters A, B and kp obtained from the fit of the thermopower 



data 



Lis] for T < 100 K by using Eq. ()34{1. In the last column we show for comparison the values 



for B' obtained from the resistivity data 



'5f\ in the range 10-100 K. 





A{^V/K^) 


kp {nm ^) 


B 


B' 


PLV film-FHNOa 


0.184±0.003 


O.lOibO.Ol 


0.132±0.002 


0.131±0.002 


HiPco fiber+H2S04 


0.083±0.007 


0.57±0.01 


0.156±0.013 


0.222±0.003 



Concerning the consistency of the fitting parameters A and kp we should make the follow- 
ing remarks. By using the values for kp shown in Table I and a simple tight binding model 
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FIG. 6: The ratio S/T as a function of temperature. The symbols are the experimental data for 
two bulks samples a, l9|. The solid lines denote the total thermopower of an individual nanotube 
obtained as explained in text. The dashed lines are the phonon-drag contributions. For clarity 
the calculated that corresponds to the PLV+HNO3 sample has been multiplied by the factor 
2. The peaks shown in the measured SjT are associated to the phonon-drag effect. 

for the estimation of the first van Hove singularity (see, for example, Ref. j49|]) the values 
we get for Ep are in good agreement with those determined from reflectivity and Raman 
measurements {q]. Also, A varies inversely with k\ in agreement with Mott's expression for 
S'^ . Moreover, the values we extract for kp support recent arguments according to which 
H2SO4 is a stronger dopant than HNO3 [9|. Namely, according to our estimation for the 



Fermi wave numbers, the Fermi level is shifted by 94 and 155 meV below the top of the 
valence band for the PLV+HNO3 and HiPco+H2S04 samples, respectively. 

In order to show clearly the effect of phonon drag in Fig. 6 we have plotted the ratio 
S/T as a, function of T. The circles and the squares are the measured values for the samples 
PLV+HNO3 and HiPco+H2S04, respectively. The dashed lines are the theoretical estimates 
for and the solid lines are the calculated values for the total thermopower. The peaks at 
T = T* are associated to phonon-drag thermopower. The shift between the experimental 
and the theoretical values for T* is due to the logarithmic term in S"^. The position of 
the peak moves towards to higher temperatures as doping increases. This dependence can 
be understood by maximizing the ratio S^/T using Eq. ( 129|) . Then we get the following 
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dependence 



T* = 1.1- 



hv.k 



k 



(36) 



B 



It is important to add that the exponential suppression of at low temperatures is 
unique for ID systems. In higher dimensions exhibits a power-law T dependence at low 



temperatures 



19 



20| . The observed peak in S/T, which is ascribed to phonon drag, underlies 



the ID character of thermopower. This adds another confirmation that thermopower in bulk 
carbon nanotube-based materials is a property of the indiyidual tube rather than a property 
of the network. 



V. CONCLUSIONS 

In summary, we have presented a rigorous model for the calculation of the phonon- 
drag thermopower in degenerately doped semiconducting SWCNTs. By using the derived 
model we investigated the dependence of on temperature, tube radius and position of 
the Fermi level. We found that decreases with the increase of the tube radius following 
approximately a law at high temperatures. In the degenerate limit, we derive a simple 
expression for which can be used as a probe for the estimation of the free carrier density 
in doped tubes. According to this expression 5*^ shows an activated T dependence at low 
temperatures. Screening effects of the carrier-phonon coupling reduce the magnitude of 
severely and result to a quasi-linear T-dependence of phonon drag at highT. Finally, we 
have compared our model with available data in acid-doped bulk samples |8|, 



9| and we found 



a very good agreement in a wide temperature range. 
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